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ABSTRACT 

The temperature and density of the hot diffuse medium pervading galaxy groups and clusters combine into 
one significant quantity, the entropy. Here we express the entropy levels and profiles in model-independent 
forms by joining two observables, the X-ray luminosity and the change in the CMB intensity due to the 
Sunyaev-Zel'dovich (SZ) effect. Thus we present both global scaling relations for the entropy levels from 
clusters and groups, and a simple expression yielding the entropy profiles in individual clusters from resolved 
X-ray surface brightness and SZ spatial distributions. We propose that our approach provides two useful tools 
for comparing large data samples with models, in order to probe the processes that govern the thermal state of 
the hot intracluster medium. The feasibility of using such a diagnostic for the entropy is quantitatively assessed, 
based on current X-ray and upcoming SZ measurements. 
Subject headings: cosmic microwave background - X rays: galaxies: clusters 



1. INTRODUCTION 

Recent X-ray observations of the hot intracluster medium 
(ICM) filling galaxy groups and clusters have indicated prob- 
lems with both the levels of the density, n ~ 10~ 3 cm -3 , and 
the distribution of the temperature T, with average values 
kT ~ 10 — 1/2 keV from very rich clusters to poor groups. 
First, the lower densities found in poor clusters and groups 
imply mass ratios m/M of ICM to dark matter (DM) con- 
siderably below the cosmic baryonic fraction approached in 
rich clusters (see Sanderson & Ponman 2003, Pratt & Arnaud 
2003). 

Second, radial profiles T(r) have been measured, if 
coarsely, in several nearby clusters; they show a roughly 
isothermal plateau extending out to substantial radii r « 0.1 - 
0.27? in terms of the virial radius R, and then a decline by 
a factor 1/2 out to R/2 (Markevitch et al. 1998, De Grandi 
& Molendi 2002, Vikhlinin et al. 2005). Although still de- 
bated in detail, this behavior differs from the predictions of 
simple modeling in terms of a polytropic equation of state 
T(r) oc « r_1 with constant index 1 < T < 5/3. It also differs 
at r < 0.2R from the outcomes of most state-of-the-art nu- 
merical simulations, where T(r) keeps rising toward the clus- 
ter center to r ~ 0.05 R as discussed, e.g., by Borgani et al. 
(2004). A central, limited dip of T(r) observed in many clus- 
ters (e.g., Piffaretti et al. 2005) and often referred to as a "cool 
core" (Molendi & Pizzolato 2001) is arguably attributed to a 
quenched "cooling flow" (see the discussion by Fabian 2004). 
Such dips involve ICM fractions of some 10~ 3 and constitute 
a specific issue not of direct concern here. 

To focus the problems in the bulk of the ICM, temperature 
and density are conveniently combined into a significant sin- 
gle quantity, the specific entropy s; here we adopt the widely 
used adiabat K = kTn~ 2 ^, which is related to s by K oc e 2s ^ k 
(Bower 1997; Balogh, Babul & Patton 1999). The quantity 
K constitutes the simplest combination of n and T that is in- 
variant under adiabatic processes in the ICM; nonadiabatic 
processes compete to produce different levels and profiles of 
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K, that cooling always lowers and various heating processes 
tend to raise. 

When the data on n and T are combined to yield K, it is 
found that (i) in moving from clusters to groups the levels 
of K decline only slowly (in fact, like K oc J 2 / 3 ) taking on 
values from around 10 3 to some 10 2 keV cm 2 ; (ii) in clusters 
the entropy profiles K(r) have to be described in terms of a 
running index 

3 a Inn 

that implies profiles K(r) flattening inwards. 

These findings may be unified under the heading of an en- 
tropy "excess" in the ICM. The baseline is provided by gravi- 
tational heating, the basic process that affects the thermal state 
of the ICM during its inflow into the structures as they are 
built up by standard hierarchical clustering (Peebles 1993). 
During the buildup, the DM component sets the gravitational 
potential wells into which the baryons fall, starting with the 
cosmic density ratio close to 0.16 (Bennett at al. 2003). If the 
baryons start cold, they fall in supersonically, and are shock- 
heated to temperatures T close to the virial value Ty (Cava- 
liere, Menci & Tozzi 1999; Tozzi & Norman 2001; Voit et 
al. 2003). However, such a process produces (see Ponman, 
Sanderson, & Finoguenov 2003): (i) entropy levels lower 
than observed in groups and poor clusters; (ii) entropy pro- 
files in rich clusters that decline too steeply inwards, in fact, 
as K(r) oc r 11 . 

Such drawbacks have raised a wide debate concerning the 
nature of the additional processes that increase the entropy 
of the ICM; suggestions include the following possibilities. 
Energy drain by radiative cooling extended well beyond the 
very central cool core may cause much of the cold gas to con- 
dense into stars, leaving a residual diffuse component with 
higher entropy (Voit & Bryan 2001); however, such a cooling 
alone would produce too much cold gas and too many stars 
(e.g., McCarthy et al. 2004). Repeated energy inputs by su- 
pernovae and active galactic nuclei preheat the gas external 
to forming groups or clusters, raise its state to a higher adi- 
abat, and are expected to hinder its inflow under any model 
(Evrard & Henry 1991; Kaiser 1991); specifically, as the in- 
flow becomes less supersonic, the shocks are weakened while 
additional entropy is carried in (Valageas & Silk 1999; Wu, 
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Fabian & Nulsen 2000; Cavaliere, Lapi & Menci 2002a). En- 
ergy impulsively discharged by powerful quasars into their 
ambient medium propagates outwards beyond the host galaxy 
into a surrounding group (Mazzotta et al. 2004a; McNamara 
et al. 2005); the propagation may occur in the form of outgo- 
ing, moderately supersonic blast waves which at the leading 
shocks raise the entropy in the gas they sweep up (Lapi, Cav- 
aliere & Menci 2005). Finally, in rich clusters thermal or tur- 
bulent diffusion may be effective in depositing and/or smear- 
ing out entropy inward of r m R (Narayan & Medvedev 2001; 
Kim & Narayan 2003; Lapi & Cavaliere, in preparation). 

Pinning down the leading process clearly requires precise 
measurements of the levels and profiles of the all-important 
adiabat K. The state variables n and T that determine K may 
be measured with X-ray observations (see Sarazin 1988). The 
bolometric bremsstrahlung emissions Lx oc n Vtr 3 w 10 42 - 
10 45 erg s" 1 primarily provides n, while spatially resolved 
spectroscopy (continuum plus emission lines) can yield the 
temperatures T when enough photons are collected. 

On the other hand, the hot ICM electrons inverse-Compton 
scatter the CMB photons crossing a cluster; this Sunyaev- 
Zel'dovich effect (SZ, Sunyaev & Zel'dovich 1972) consti- 
tutes a major ICM (and cosmological) probe in the ^wave 
and submm bands, as reviewed by Rephaeli (1995), Birkin- 
shaw (1999), and Carlstrom, Holder & Reese (2002). The 
intensity of the thermal SZ effect is given by the Comptoniza- 
tion parameter y oc nT R, which constitutes a good measure 
of the electron pressure p oc n T in the ICM; so far it has been 
observed in many rich clusters at the expected levels y~ 10" 4 
(Zhang & Wu 2000; Grainge et al. 2002; Reese et al. 2002; 
Benson et al. 2004). The main dependencies of L x oc n 2 and 
y oc p strongly motivate us to employ the results of X-ray and 
SZ observations combined as to provide the adiabat in the 
form K oc prT 5 ! 3 , in order to better probe properties of, and 
processes in the ICM. 

2. SCALING THE ENTROPY LEVELS FROM CLUSTERS TO 
GROUPS 

The levels Ko.i of the entropy (usually sampled at r w 
0. 1 #200 ~ 0.08 R, see Ponman et al. 2003) and of the X-ray lu- 
minosity L x (dominated by the inner ICM anyway) are linked 
by the relation 

K 0A /K g oc (L x /L g r^ 3 (T/T v ) 7 / 6 , (2) 

derived in Appendix A. In moving from clusters to groups, 
this highlights the deviations from the quantities K g oc 
H~ 4 / 3 Ty and L g oc H Ty that provide the baseline scaling - 
including the Hubble parameter H(z) - set by the pure gravita- 
tional heating (see § 1); the latter implies T w Ty, and baryon 
densities simply proportional to those of the DM (Kaiser 
1986). 

Eq. (2) has been first proposed by Cavaliere, Lapi & Menci 
(2002b); similar expressions have been subsequently used by 
Dos Santos & Dore (2002), and McCarthy et al. (2003), in 
the context of specific models for entropy enhancement. Here 
we stress the model-independent nature of Eq. (2); the latter 
holds whenever the ICM fills in a statistical, virial-type equi- 
librium a DM potential well, and depends weakly on specific 
DM distributions or specific equations of state for the ICM. 
This is demonstrated in Appendix A and Table Al, where we 
show that from clusters to groups the prefactor on the r.h.s. 
of Eq. (2) depends weakly on T, so as to closely preserve the 
above scaling. 



We illustrate the role that Eq. (2) may play in probing the 
processes that set the ICM thermal state. We first recall that 
in poor clusters and groups with kT ranging from about 4 to 
about 1/2 keV the observed values of L x are found to be pro- 
gressively lower than L g by factors from 10 _1 to 10~ 3 , with 
a steep decline close to Lx oc T 3 (Osmond & Ponman 2004). 
In other words, moving toward smaller and cooler structures 
the ICM is found to be progressively underluminous, hence 
underdense relative to the gravitational values. The out- 
come is even more surprising on two accounts: smaller struc- 
tures ought to have condensed earlier from a denser Universe; 
moreover, emission lines from highly excited metals add to 
the bremsstrahlung continuum to yield a flatter gravitational 
scaling L x w L g oc T for kT < 2 keV (see Appendix A). Ac- 
cording to Eq. (2), the systematically lower values of Lx /L g 
will correspond to higher entropy levels, exceeding K g by fac- 
tors up to 10. In detail, the steep relation L x oc T 3 will reflect 
into a flat behavior Kq i oc T 2 ! 3 , which indeed agrees with re- 
cent data analyses (see Ponman et al. 2003). 

On the other hand, a number of authors (e.g., Mahdavi 
et al. 1997; Roussel, Sadat & Blanchard 2000) have ques- 
tioned the significance of the data concerning groups, given 
various observational and systematic uncertainties related to 
the limited statistics and to the necessary subtraction of con- 
siderable contributions from single galaxies. Recent critical 
assessments of the data (Mushotzky 2004; Osmond & Pon- 
man 2004) acknowledge the relevance of these problems in 
individual groups, but confirm the general trend toward lower 
densities in poorer systems, albeit with a wide scatter. 

To cross-check these low X-ray luminosities Lx, one may 
preliminarily use the SZ effect; in fact, these two independent 
probes of the ICM density are expected to correlate according 
to 

y /y g cx (L x /L g ) l l\T /T v ) 3 / 4 . (3) 
Here yo is the Comptonization parameter integrated along a 

3 /2 

"central" line of sight (l.o.s.), and y g oc HT V ' gives its grav- 
itational scaling after Cole & Kaiser (1988); see Appendix 
A for the derivation and added detail, including a discus- 
sion of the contributions to yo from the outer ICM. When 
L x oc r 3 applies, Eq. (3) predicts the scaling yo oc T 2 to hold 
in equilibrium conditions (Cavaliere & Menci 2001). Such 
systematically lower values of yo/y g oc T l l 2 may be tested or 
bounded with high-sensitivity SZ observations of poor clus- 
ters and groups, as initiated by Benson et al. (2004). 

Then one can proceed to directly obtain the entropy levels 
from the equilibrium relations Eqs. (2) and (3). Upcoming SZ 
measurements will provide high-sensitivity data, which will 
be conveniently joined with X-ray fluxes to obtain entropy 
levels that scale as 

Ko.rocH-Wy^Lx 2 , (4) 
in equilibrium conditions but otherwise in a model- 
independent way. In moving to groups where line emission 

is important, Eq. (4) goes over to Ko.i oc H~ 5 I 3 y 4 / 3 ' L x l ; Ap- 
pendix A gives the details, including the slow dependence on 
T of the prefactor. Note that the present approach bypasses 
two difficulties (discussed, e.g., by Mazzotta et al. 2004b) in 
deriving the ICM temperatures from spatially resolved X-ray 
spectroscopy: (i) the latter requires more photons than imag- 
ing, so with comparable observation times and efficiencies, 
precise T are obtained at smaller radii than n; ii) in a multi- 
phase medium spectroscopic determinations of T are biased 
toward the colder regions. 
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The integrated SZ flux Y oc y^R 2 jd\, where <aU(z) is the an- 
gular diameter distance, may be used instead of the central 
parameter yo\ in terms of Y, Eq. (4) may be recast as 



Ko.^df^H^Y^L^ 9 



(5) 



Using this relation has two advantages over Eq. (4). Not only 
is Y a more robust observable in clusters with a cool core (e.g., 
Benson et al. 2004), but also Y oc nT R 3 has the same R 3 de- 
pendence as L x . So the combination is largely free of the size 
effects due to the scaling R oc T l l 2 for coeval structures (cf. 
Eq. [A2]); it rather highlights the dependence on the baryonic 
fraction. 

So far we have dealt with the overall entropy levels deter- 
mined from the integrated X-ray and SZ observables. But the 
entropy profiles of individual clusters carry additional infor- 
mation; clearly this will be difficult to extract from the prefac- 
tors in global relations like Eqs. (2)-(5), encased as it is under 
integrated forms within ratios or products of shape factors, as 
made explicit in Appendix A. 

3. RESOLVING THE ENTROPY PROFILES IN CLUSTERS 

To go beyond the above global scaling laws, spatially re- 
solved X-ray and SZ observations are clearly required. In X 
rays, Chandra and XMM-Newton are already providing data 
at a few arcsec resolutions; but obtaining both the good spa- 
tial and spectral resolution needed to measure T(r) and derive 
K(r) still presents a challenge, especially in the outer ICM of 
distant clusters. Improvements will again be provided by join- 
ing measurements of X-ray brightness and of SZ effect with 
comparable levels of spatial resolutions. 

In fact, the prospective good resolution and high sensitiv- 
ity (added to multi-frequency capabilities) of the many up- 
coming SZ experiments (for a recent review, see Birkinshaw 
2004) stimulate renewed impetus for considering the depro- 
jected quantities that have been discussed from time to time 
(e.g., Silk & White 1978, Cavaliere 1980). 

Here our specific proposal is to use joint X-ray and SZ mea- 
surements to derive directly from the data the radial entropy 
profiles K(r) in the ICM. We will work in terms of the X-ray 
surface brightness and of SZ intensity distributions, integrated 
along a l.o.s at a projected distance s from the cluster center. 
For spherically averaged ICM distributions these observables 
are related to the corresponding volume quantities by 



W= 9 nl ^4 f dr^=n\r)^f{7) 
27r(l+z) 4 J s Vr 2 -s 2 



2 koj f r 
y(s)= 2 dr—==n(r)T{r) , 



(6) 



where the constants are defined after Eq. (Al). Note how 
strongly the X-ray surface brightness decreases as the redshift 
z increases, while the SZ effect is z-independent, apart from 
the dilution effect related to the size of the instrumental beam. 

To construct the adiabat K = kT n" 2 / 3 we need to retrieve the 
volume quantities n(r) and T(r) from the projected quantities 
£ x (s) and y(s) given by Eqs. (6). Actually, the latter are in the 
form of two Abel integral equations, which may be solved for 
the volume quantities along the lines recalled in Appendix B; 
the results come to 



n 2 Vkf-- 



nkT- 



4£ 1 / 2 (l+z) 4 1 d [ r J s 



A 



\_d_ r 

r dr J R 



ds 



Vs 2 ^? 



£ x (s) = £ x (r) 



(7) 



m e c 2 1 d f r 



-— ds 



-k/jt i" dr J R - r 2 



y(s) = y(r) 



We stress that the two volume quantities £ x (r) and y(r) are 
defined directly in terms of the observed surface or l.o.s. ob- 
servables £ x (s) or y(s). 

Eqs. (7) are easily combined (over the spatial range com- 
mon to both data sets) to yield the entropy profile in the sim- 
ple, model-independent form 



K(r) = y 14 / 9 (r)£ x W/9 (r) 



(8) 



So far we have not assumed any specific equation of state for 
the ICM, rather we derive it from K(r) above; in particular, 
the local polytropic index defined in Eq. (1) reads 



rw = 3 



i dln£ x 
' d\ny 



1 



(9) 



To illustrate the importance of the profile K(r) we contrast 
two physical conditions. On the one hand, a reference base- 
line is provided once again by the gravitational heating asso- 
ciated with smooth accretion of successive shells (involving 
DM plus gas, as started by Bertschinger 1985) during the hi- 
erarchical buildup of a rich cluster. This "surface" process 
taking place at the outskirts of the growing structure yields 
K(r) oc r l 1 and T w 1.1 in the region r R (Tozzi & Nor- 
man 2001). On the other hand, the actual inflow is likely 
to include additional "volume" processes induced by infall 
and sinking of lumps, the well-known merging component in 
the process of structure formation (see Norman 2005). These 
dynamical events can drive subsonic turbulent diffusion and 
related entropy transport from the outskirts toward the inner 
regions (see Goldman & Rephaeli 1991; Ricker & Sarazin 
2001; Inogamov & Sunyaev 2003); clearly these processes 
will enforce entropy profiles flatter than K(r) oc r 11 , a matter 
expanded upon elsewhere (Lapi & Cavaliere, in preparation). 
So it is seen that the entropy profiles can help unveiling the 
actual physical conditions in clusters. 

Note that Eqs. (7) also yield the density and temperature 
profiles in the form 



n{r) = cf{r)y- x l\r) 



, x y,,j y,, kT(r) = y A l\r)C x '\r) . (10) 

Correspondingly, the local cooling timescale for 
bremsstrahlung emission reads 



tc(r) = 



3k l l 2 y(r) 
A £x(r) 



(11) 



Comparing this to the Hubble time provides a criterion to 
avoid the inclusion of a region like a cool core, where the 
physical conditions are very different from the rest of the clus- 
ter. 

So far, our derivations have been free of specific assump- 
tions concerning models for the DM gravitational potentials, 
or hydrostatic equilibrium, or equation of state. Assuming 
now detailed hydrostatic equilibrium, one can use the ICM 
quantities to trace the DM mass distribution (e.g., Bahcall & 
Sarazin 1977; Fabricant & Gorenstein 1983), and obtain the 
result 

M ^ r 2 d y4/3^ y . 

r ~ AG^m p £ 2 / 3 (r) dr 
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FIG. 1 . — Input profiles for the mock experiments described in § 4: X-ray surface brightness, SZ intensity distribution, and radial profiles of density, temperature, 
entropy and polytropic index. Dot-dashed lines are for the simple polytropic model A with T = 1.1, while solid lines are for the nonpolytropic model B with 
running index T(r); see main text for details. 



here /i w 0.6 is the mean molecular weight for a plasma of 
one-third solar metallicity, m p is the proton mass and G is 
the gravitational constant. Equivalently, Eq. (12) may be used 
to derive a model for the DM gravitational potential, to be 
compared with specific forms such as King's (1972) and NFW 
(Navarro, Frenk & White 1997), or with the data obtained 
from gravitational lensing (e.g., Lombardi et al. 2005). 

4. MOCK EXPERIMENTS 

Given that the imaging capabilities of the present X-ray in- 
struments (Chandra and XMM-Newtori) have attained resolu- 
tions down to about 1" and fractional sensitivities of order 
10" 3 or better, what can be learned right now by joining these 
X-ray data with the current, considerably coarser SZ observa- 
tions? 

In fact, here problems arise even for the global scaling laws, 
and more for profile reconstructions. As discussed in detail 
by Diaferio et al. (2005), considerable discrepancies as to the 
scaling laws are found among current SZ samples produced 
with different instruments, even when the samples have many 
sources in common. On the other hand, to effectively con- 
strain a fit to the observed SZ distributions often it has been 



found necessary to use also X-ray information (see Reese et 
al. 2003, Benson et al. 2004); this breaks the independence 
of these two observables that is important for deriving K(r). 

A step forward has recently been taken by the observations 
of the X-ray cluster RX J 1347-1 145 with the Nobeyama tele- 
scope (Komatsu et al. 2001); the good resolution of order 
15" obtained in these measurements enabled deriving a truly 
independent SZ profile. In analyzing these data Kitayama et 
al. (2004) excise the most asymmetric sector, and spherically 
average the rest to deconvolve the temperature profile, find- 
ing results consistent with - but still no better than - the X-ray 
spectroscopy. 

In the light of upcoming SZ capabilities of sub-arcmin res- 
olution and sensitivities better than 10 /iK, we consider next 
the prospects from measurements with instruments or experi- 
ments such as SZA (Holder et al. 2000), AMI (Jones 2002), 
AMiBA (Lo 2002), MITO (Lamagna et al. 2005), OLIMPO 
(Masi et al. 2003), ACT (Kosowsky 2003), SPT (Carlstrom 
et al. 2003), APEX-SZ (Schwan et al. 2003), CARMA 
(Woody et al. 2004). It is expected that many hundreds 
of clusters and groups will be detected in a survey mode, 
with pointed observations of many clusters to sensitivity lev- 
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FIG. 2. — Reconstructed profiles from the the mock experiments of § 4: radial profiles of density, temperature, entropy, and polytropic index are shown from 
top to bottom as solid lines, with \a error bars; dashed lines show the input profiles. Left column refers to model A and right column to model B. 



els of 10 fiK or better. In preselected areas, resolutions down 
to a few arcsecs will be eventually provided by ALMA (see 
|http : / /www . alma . nrao . edu /), at frequencies on both 
sides of 220 GHz, the crossover point where the thermal SZ 
effect goes from negative to positive. 

We build mock renditions of the prospective observations 
as described next. We make use of two independent, flexi- 
ble parameterizations for the spatial distributions of the X-ray 
brightness and SZ effect; both are of the form 



/(s) = /(0) 



1 + 



(13) 



often adopted for either one, though best suited to clusters 
with no cool core. Here a is the inner, characteristic radius 
and b is the shape parameter for the outer profile (e.g., Reese 



et al. 2002; Schmidt, Allen & Fabian 2004). Note that a/Vb 
provides an effective "core radius", while I(s) oc s~ 2b approx- 
imately holds for large r. To preserve the independence of 
the two fits, we allow for different values of the parameters 
o-XiO-sz an d bx,bsz in Eq. (13); note that the isothermal (3- 
model (Cavaliere & Fusco-Femiano 1976; Jones & Forman 
1984) would obtain as the particular case where ax = asz, 
b x = 3/3- 1/2 re 1.5 and fc sz = (3 ( 3-l)/2«0.5 apply. 

First, let us disregard the limitations due to finite reso- 
lution and precision; then the fitting formulae Eq. (13) are 
easily Abel-inverted, and lead to sharp volume distributions 
n 2 T l l 2 oc [l+(r/a x ) 2 T hx - 1/2 and nT oc [1 +(r/a sz ) 2 r isz ~ 1/2 
for the X-ray emissivity and SZ effect, respectively. Using 
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these in Eqs. (8) and (9), the entropy profile reads 



5(l+2fe)/9 



K(r)/K(0)=[l + {r/a x ) 
while the polytropic index is given by 



\ + (r/a sz )' 



-7(l+2b sz )/9 



(14) 



r(r) = 3 



2 + 4b x r 2 + a 2 s 



sz 



\+2b sz r 2 + a 2 



1 



(15) 



These runs are plotted in Fig. 1 for the specific parameter val- 
ues given in the caption. Note that, when a x = «sz is set, a 
constant value of T = (3 + 6bsz)/(4-b x -2bsz + 1) would ob- 
tain. Then the bounds (2b x - l)/4 < b sz <(5b x -\)/l ap- 
ply; the upper limit corresponds to T = 5 /3, i.e., the isentropic 
profile, while the lower limit corresponds to T = 1, i.e., the 
isothermal profile. 

Next we proceed to generate and analyze mock SZ experi- 
ments on introducing prospective resolutions and precisions. 
We focus on two input models: A) a simple polytropic model 
with constant index T s=a 1.1, which corresponds to Eq. (13) 
with bx = 2.1, bsz = 0.9, a x j\fb x = 0.08/? and a x = asz', 
B) a nonpolytropic model based on Eq. (13) with b x = 2.1, 
bsz = 0.9, a x / ^/b x = 0.\5R and a S z/Vbsz = 0.21 R, implying 
a running index T(r). These input profiles are illustrated in 
Fig. 1. 

Having in mind rich clusters at z < 0.1, we bin the data 
into logarithmically equal intervals limited by a central reso- 
lution of 10"; higher resolution, such as eventually attainable 
with ALMA, will extend the limiting z beyond 0.5. Follow- 
ing Yoshikawa & Suto (1999; see also Tsutsumi et al. 2005, in 
preparation), we assign 1 a uncertainties Al x w 2 x lO^I^O) 
and Alsz ~ 10~ 2 /sz(0) to the input X-ray surface brightness 
and SZ intensity distributions, respectively. Using a stan- 
dard Monte Carlo scheme we randomly sample Gaussian- 
distributed values for I x and I$z, then perform the reconstruc- 
tion, and finally average over 1000 realizations. 

Fig. 2 shows for each model A and B the input {dashed 
lines) and reconstructed profiles (solid lines with 1 a error 
bars) of the ICM density n(r), temperature T(r), entropy K(r) 
and polytropic index T(r). 

It is seen that the average reconstructed profiles are quite 
similar to the input profiles over a considerable range, in spite 
of uncertainties and of some systematic deviations increas- 
ingly induced at small and large radii by our emulated reso- 
lution and sensitivity. It is also seen that the uncertainties in 
T(r), dominated by the assumed limitations of the SZ obser- 
vations, considerably exceed those in n(r), mainly contributed 
by the X-ray observations. On the other hand, the former do 
not grow fast toward the dimmer outskirts, unlike the X-ray 
observations particularly of distant clusters in the presence of 
instrumental or astrophysical background. We add that un- 
certainties and deviations reduce to about one half when we 
assume levels of sensitivity and of spatial resolution twice 
higher than the values used in Fig. 2. 

5. DISCUSSION AND CONCLUSIONS 

To make optimal use of the precise and resolved, upcom- 
ing SZ data we have proposed a simple formalism, which is 



model-independent even though formulated here for spher- 
ically averaged symmetry. Different geometries may be 
considered along the lines discussed by Binney & Strimpel 
(1978), and advanced by Zaroubi et al. (2001). 

We propose that SZ and X-ray measurements jointly pro- 
vide useful diagnostic tools aimed at assessing the thermal 
state of the ICM with the use of large data sets. Specifically, 
Eqs. (4) or (5) provide simple, global scaling laws for the lev- 
els of the central entropy Km from clusters to groups; on the 
other hand, in relatively nearby rich clusters Eqs. (8) and (9) 
directly yield the resolved profiles K(r) and the running poly- 
tropic index T(r). The former relations will provide insight 
on the processes governing the overall thermal state of the 
ICM in clusters and in groups. The latter relations will probe, 
e.g., the effectiveness of thermal or turbulent diffusion in es- 
tablishing the local equation of state for the ICM within rich 
clusters. 

We stress that the scaling in Eq. (3) yields what actually rep- 
resents a lower limit to yo, valid for the ICM in equilibrium 
within the potential wells. As shown in Table Al and com- 
mented below it in Appendix A, the contribution to yo from 
the outer range r > 0.2R is limited to 10- 15%, and encased 
in the shape factors; their slow dependence on system tem- 
perature can only yield moderate deviations when the outer 
and inner regions are described by the same monotonically 
decreasing hydrostatic distribution. By the same token, up- 
ward deviations of yo will strongly indicate contributions from 
outer nonequilibrium plasma, overpressured and outflowing 
from the wells, an issue that discussions with P. Mazzotta 
stimulated us to stress. In fact, Lapi, Cavaliere & De Zotti 
(2003, see their Fig. 3) have computed the transient upper 
limits expected in 5- 10% early bright galaxies and surround- 
ing groups when powerful quasars flare up and drive outgoing 
blast waves that sweep out part of the ambient medium while 
raising its pressure and entropy. 

In summary, SZ observations with sensitivities of 10^K or 
better and sub-arcmin resolutions will provide the means to 
scale entropy levels in poor clusters and groups out to z ~ 1 .5. 
Resolutions of order 10" will yield pressure and entropy pro- 
files, most valuable in the outskirts of relatively nearby clus- 
ters as long advocated by Rephaeli (e.g., 2005). Higher reso- 
lutions of a few arcsecs will enable reconstructing the entropy 
profiles in rich clusters toward z ~ 1, where SZ may effec- 
tively compete with X rays being unaffected by cosmological 
dimming. These SZ effect frontiers will afford unprecedented 
insight on the astrophysics of the ICM. 
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APPENDIX 
A. THE SCALING LAWS DERIVED 
The integrated X-ray luminosity and the "central" SZ parameter (in the non-relativistic case) are given by 
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m e c z J n 2 T 2 



where we have used averaged spherical symmetry in terms of the non-dimensional radial coordinate x = r/R, and have normalized 
the electron density n and the temperature T to their values n 2 and T 2 at r = R. Here the expression Ao T e approximates the X-ray 
emissivity in terms of the normalization Ao (taking on values around 2 x 10~ 27 erg cm 3 s" 1 K" 1 / 2 for bremsstrahlung emission) 
and of the power-law index e (ranging along the cooling curve from e w 1/2 for bremsstrahlung-dominated emission to e w — 1/2 
for important line emissions at kT < 2 keV; see Sutherland & Dopita 1993). 

The statistical equilibrium of the ICM within a DM potential well of depth marked by the virial temperature Ty implies the size 
dependence 

T v <xH 2 R 2 . (A2) 

Here H(z) = [0^(1 + z) 3 + QaY^ 2 is defined in terms of the virialization redshift z, and of the present density Q M ~ 0.27 for the 
non-relativistic matter (DM + baryons) and Q\ w 0.73 for the dark energy in the Concordance Cosmology (Bennett et al. 2003); 
the approximation H(z) cx (1 +z) holds for z < 1. 
Thus we obtain the scaling laws 

L x cx S X H~ 3 {T v /T 2 fl 2 n 2 T 2 e+3/2 y cx S y H~ x {T v /T 2 ) l l 2 n 2 T 2 3/2 , (A3) 

in terms of the two shape factors S x , S y defined by the integrals appearing in Eqs. (Al). As to the adiabat, its levels A"o.i at 
/■kO.IR may be written as 

in terms of the other shape factor Sk- We will see below that all factors: T 2 /Ty, S x , S y , and Sk, vary only weakly from clusters 
to groups, and so constitute minor corrections to the explicit scaling laws above. This stems from the circumstance that the 
structural function n(r)/n 2 entering the shape factors depends weakly on T 2 from clusters to groups, as observed by Sanderson & 
Ponman (2003) and Pratt & Arnaud (2003). Weak dependence of n(r) is also expected from all equilibrium models (as discussed 
in depth by Cavaliere et al. 2002b); in fact, these show a weaker and weaker dependence in passing from King (1972) to NFW 
gravitational potentials (due to their DM concentration rising from clusters to groups, see Navarro, Frenk, & White 1997), or 
from isothermal to polytropic ICM (due to a lesser density rise toward the center). 
On eliminating n 2 from Eqs. (A3) and (A4) we obtain the relations 

yo/y g cx (T v /T 2 y 1+ ^ 2 (L x /L g ) 1/2 K 0A /K g cx (T v /T 2 y 1 -^ 3 (L x /L g )~ l/i , (A5) 

which for e = 1 /2 yield Eqs. (2) and (3) of the main text. We have normalized Ko. i , yo andL x to the respective gravitational scaling 

K g cx H~ A I 3 T v , y g cx HTy 2 , L g cx HTy +3 ^ 2 . Recall from § 2 that these obtain on considering T 2 w T v , and n 2 cx H 2 following the 
DM density, in turn proportional to the background's; note that when line emission is important e = -1/2 and L g cx T v apply, as 
stated in § 2. 

Moreover, on eliminating n 2 and the main T 2 dependence from Eqs. (A3) we find 

Ko.! OC // -(18-4 e )/(9-6 e ) }; (18 + 4 e )/(9-6 e ) L -4/(3-2 e ) ^ (A6) 

which for e = 1/2 yields Eq. (4) of the main text; for simplicity we have omitted the slowly T-dependent prefactors. In full, these 

read (T v /T 2 T 1+e / 2 S y S x 1/2 , (T v /T 2 )~ l ~ e / 3 S K S x 3 and (T v / T 2 )®~ 2 e) /< 9 - 6 «) S x ^ 3 ~ 2 e) <Sy (18+4 e) / (9_6 e) S K on the r.h.s. of Eqs. (A5) and 
(A6), respectively. 

For the purpose of illustrating the slow dependence of Ty/T 2 , S x , S y , and Sk we now assume a specific DM potential well, and 
detailed hydrostatic equilibrium of the ICM with a specific relation between n(r) and T(r). For the latter we adopt a polytropic 
distribution with T(r) cx n(r) r_1 in terms of an average, constant index r w 1.1 — 1.3; then the temperature run writes T(r)/T 2 = 
l + (r— 1) f3 A<f)(r)/T, in terms of the potential difference A<fi(r) = [<£>(/?)- $(r)]/er 2 normalized to the DM one-dimensional 
velocity dispersion a. As to the DM potential $(r), we use the simple King form $ = -9er 2 ln[(r/r c ) + (l +r 2 /r 2 )'' 2 ] r c /r with 
r c = R/12, but similar or better results obtain on using the NFW model. The parameter (3 = [im p a 2 /kT 2 = Ty/T 2 encodes the T 2 
dependence of the normalized profiles; it takes on values varying in the range from 0.7 in rich clusters to 0.5 in poor groups. 

In Table Al we show our results; for simplicity e = 1/2 is adopted, but we add that S x changes less than 15% for e = -1/2 (e.g., 
for T = 1.2 the value decreases from 0.73 to 0.65 in groups with f3 = 0.5). To put these changes in context, recall from § 2 that 
in moving from rich clusters with kT w 10 keV to groups with kT m 1 keV, L x cx T 3 is observed to decrease by factors 10 2 or 
more, Ko.i cx T 2 I 3 by at least a factor 5, and yo cx T 2 is expected to decrease by 30 or more; thus the T dependencies of the shape 
factors in the scaling relations constitute minor corrections. Actually, even slower variations occur in the prefactors of Eqs. (A5) 
and (A6) as they contain combinations of shape factors whose variations tend to compensate; e.g., for T = 1.2 the prefactors in 
Eqs. (A5) vary by 0.9 and 1.5, respectively, from clusters to groups. 
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TABLE Al 
Shape factors, dependence on T 2 







r = i.i 






r = i.2 






r = 1.3 







S x 


Sy 


S K 


Sx 


Sy 


S K 


Sx 


Sy 


Sk 


0.7 


1.66 


7.43 


0.19 


1.17 


5.78 


0.33 


0.96 


4.94 


0.47 


0.6 


1.15 


5.29 


0.24 


0.91 


4.42 


0.37 


0.79 


3.93 


0.51 


0.5 


0.84 


3.82 


0.29 


0.73 


3.37 


0.43 


0.66 


3.10 


0.56 



NOTE. — The shape factors are computed for a polytropic ICM in equilibrium within a King's DM potential well. Recall that /3 = Ty/T 2 ranges from 0.7 to 
0.5 in moving from rich clusters to poor groups. 



As to sensitivity of the SZ effect to the outer wings of the ICM distribution, note the bounds to their contribution set by the 
observed outer decline of the temperature T(r); to give an example, the contribution to yo from the outer range r > Q.2R comes to 
15 - 10% when the polytropic index ranges between V w 1 .1 - 1.3. The T-dependence from clusters to groups, the matter relevant 
to the scaling relations Eqs. (3)-(5), is bounded by the slow variation of the overall shape factor S y given in Table Al. 

B. THE ABEL EQUATION 
We recall that the Abel integral equation has the form 

f Fit) 

f(x)= / dt —r== ; (Bl) 

here f(x) is a known function, q a constant (that may be large), and F(t ) the unknown function. To derive the latter, multiply both 
sides of Eq. (Bl) by (x-£)~ 1//2 and integrate overx in the interval [£, q\. 

where the second equality is obtained on exchanging the integration order; the last integral is simply tt. Finally, differentiate both 
sides of Eq. (B2) with respect to £ to obtain 

rS fix) 



F ^)=-4- c I dx-L^=, (B3) 

71" d£ J q 



which is the Abel inversion formula. 

Note that on differentiating with respect to £ (with due attention paid to the singularity of the kernel), Eq. (B3) may be recast 
into the alternative form 

i™ 1 f*s f ' (x) JML ™ 
F(0 = - dx + - , (B4) 

T Jq Vx-£, 7T VQ-Z 

where the prime denotes differentiation with respect to x. But, apart from problems with the boundary term f(q) /tt 
Eq. (B4) is less useful in the context of the main text since it involves differentiation of the data or of their fit, with the related 
uncertainties (possibly enhanced by real dumpiness in the data) as discussed in detail by Lucy (1974) and Yoshikawa & Suto 
(1999). 
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